MITTAG-LEFFLER FUNCTIONS AND THEIR APPLICATIONS 



H.J. HAUBOLD 
Office for Outer Space Affairs, United Nations 
Vienna International Centre, P.O. Box 500, A-1400 Vienna, Austria, 
and Centre for Mathematical Sciences, Pala Campus, 
Arunapuram P.O., Pala, Kerala-686574, India. 
hans.haubold@unvienna.org 

A.M. MATHA1 

Department of Mathematics and Statistics, McGill University 
Montreal, Canada H3A 2K6 
and Centre for Mathematical Sciences, Pala Campus, 
Arunapuram P.O., Pala, Kerala-686574, India. 
mathai@math . mcgill . ca 

R.K. SAXENA 

Department of Mathematics and Statistics, Jai Narain Vyas University 
Jodhpur-342005, India. 
ram.saxena@yahoo.com 

Abstract 



Motivated essentially by the success of the applications of the Mittag-Leffler functions in many areas of 
science and engineering, the authors present in a unified manner, a detailed account or rather a brief survey 
of the Mittag-Leffler function, generalized Mittag-Leffler functions, Mittag-Leffler type functions, and their 
interesting and useful properties. Applications of Mittag-Leffler functions in certain areas of physical and 
applied sciences are also demonstrated. During the last two decades this function has come into prominence 
after about nine decades of its discovery by a Swedish Mathematician CM. Mittag-Leffler, due its vast 
potential of its applications in solving the problems of physical, biological, engineering and earth sciences etc. 
In this survey paper, nearly all types of Mittag-Leffler type functions existing in the literature are presented. 
An attempt is made to present nearly an exhaustive list of references concerning the Mittag-Leffler functions 
to make the reader familiar with the present trend of research in Mittag-Leffler type functions and their 
applications. 

1. Introduction 



The special function 

00 k 

and its general form 

00 p. 

E a ,p(z) = wTTT 1 \ 1 a i R(a) > 0, R(/?) > 0, z e C (1.2) 

with C being the set of complex numbers are called Mittag-Leffler functions (Erdelyi, et al., 1955, Section 
18.1). The former was introduced by Mittag-Leffler (Mittag-Leffler, 1903), in connection with his method 
of summation of some divergent series. In his papers (Mittag-Leffler, 1903, 1905), he investigated certain 
properties of this function. The function defined by (1.2) first appeared in the work of Wiman (Wiman, 
1905). The function (1.2) is studied, among others, by Wiman (1905), Agarwal (1953), Humbert (1953) 
and Humbert and Agarwal (1953) and others. The main properties of these functions are given in the book 
by Erdelyi et al. (1955, Section 18.1) and a more comprehensive and a detailed account of Mittag-Leffler 
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functions are presented in Dzhcrbashyan (1966, Chapter 2). In particular, the functions (1.1) and (1.2) are 
entire functions of order p = 1/a and type a = 1; see, for example, (Erdelyi, et al., 1955, p. 118). 

The Mittag-Lcfflcr function arises naturally in the solution of fractional order integral equations or 
fractional order differential equations, and especially in the investigations of the fractional generalization 
of the kinetic equation, random walks, Levy flights, super-diffusive transport and in the study of complex 
systems. The ordinary and generalized Mittag-Leffler functions interpolate between a purely exponential 
law and power-law like behavior of phenomena governed by ordinary kinetic equations and their fractional 
counterparts, see Lang (1999a, 1999b), Hilfer (2000), Saxena et al. (2002). 

The Mittag-LefHer function is not given in the tables of Laplace transforms, where it naturally occurs 
in the derivation of the inverse Laplace transform of the functions of the type p a (a + bp 13 ), where p is the 
Laplace transform parameter and a and b are constants. This function also occurs in the solution of certain 
boundary value problems involving fractional integro-differcntial equations of Volterra type (Samko et al., 
1993). During the various developments of fractional calculus in the last four decades this function has 
gained importance and popularity on account of its vast applications in the fields of science and engineering. 
Hillc and Tamarkin (1930) have presented a solution of the Abel- Volterra type equation in terms of Mittag- 
LefHer function. During the last 15 years the interest in Mittag-LefHer function and Mittag-LefHer type 
functions is considerably increased among engineers and scientists due to their vast potential of applications 
in several applied problems, such as fluid flow, rheology, diffusive transport akin to diffusion, electric networks, 
probability, statistical distribution theory etc. For a detailed account of various properties, generalizations, 
and application of this function, the reader may refer to earlier important works of Blair (1974), Bagley and 
Torvik (1984), Caputo and Mainardi (1971), Dzherbashyan (1966), Gorenflo and Vessella (1991), Gorenflo 
and Rutman (1994), Kilbas and Saigo (1995), Gorenflo et al. (1997), Gorenflo and Mainardi (1994, 1996, 
1997), Gorenflo, Luchko and Rogosin (1997), Gorenflo, Kilbas and Rogosin (1998), Luchko (1999), Luchko 
and Srivastava (1995), Kilbas, Saigo and Saxena (2002, 2004), Saxena and Saigo (2005), Kiryakova (2008a, 
2008b), Saxena, Kalla and Kiryakova (2003), Saxena, Mathai and Haubold (2002, 2004, 2004a, 2004b, 2006), 
Saxena and Kalla (2008), Mathai, Saxena and Haubold (2006), Haubold and Mathai (2000), Haubold, Mathai 
and Saxena (2007), Srivastava and Saxena (2001), and others. 

This paper is organized as follows: Section 2 deals with special cases of E a (z). Functional relations of 
Mittag-LefHer functions are presented in Section 3. Section 4 gives the basic properties. Section 5 is devoted 
to the derivation of recurrence relations for Mittag-Lcfflcr functions. In Section 6, asymptotic expansions 
of the Mittag-Leffler functions are given. Integral representations of Mittag-Lcfflcr functions are given in 
Section 7. Section 8 deals with the H-function and its special cases. The Mellin-Barnes integrals for the 
Mittag-LefHer functions are established in Section 9. Relations of Mittag-Leffler functions with Riemann- 
Liouville fractional calculus operators are derived in Section 10. Generalized Mittag-Leffler functions and 
some of their properties are given in Section 11. Laplace transform, Fourier transform, and fractional 
integrals and derivatives are discussed in Section 12. Section 13 is devoted to the application of Mittag- 
Leffler function in fractional kinetic equations. In Section 14, time-fractional diffusion equation is solved. 
Solution of space- fractional diffusion equation is discussed in Section 15. In Section 16, solution of a fractional 
reaction-diffusion equation is investigated in terms of the H-function. Section 17 is devoted to the application 
of generalized Mittag-Leffler functions in nonlinear waves. Recent generalizations of Mittag-LefHer functions 
are discussed in Section 18. 

2. Some special cases 

We begin our study by giving the special cases of the Mittag-LefHer function E a (z). 



E (z) = 



1 



1 - z 



(ii) 



E 1 (z)=c z 



(iii) 



E 2 {z) = cosh(Vz),z e C 
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(iv) 



E 2 (-z 2 ) = cosz, z e C 



(v) 



e* s + 2e" 



cos 



V3 



-Z3 



,zeC 



(vi) 



Ei{z) = \ cos(zi ) + cosIi(zi ) 



,ze C 



(wm) £i(±z2) =c z l + crf(±z2) = c z erfc(=Fz2),z e C 

where erfc denotes the complimentary error function and the error function is defines as 

2 f z 

erf (z) = —= / exp(-i 2 )dt, erfc(z) = 1 - erf (z), z e C. 
V 71 " Jo 

For half- integer n/2 the function can be written explicitly as 



(viii) 



E*(z) = F n _ 1 



12 n - 1 z 2 



+ - 



2 (n+l)/2 z 



l^n-i 1; 



n + 2 n + 3 in z 2 



(ix) 



E\,1\Z) = , hj2,7\Z) = = . 



3. Functional relations for the Mittag-Leffler functions 

In this section, we discuss the Mittag-Leffler functions of rational order a — m/n, with m,n € N 
relatively prime. The differential and other properties of these functions are described in Erdelyi, et al.(1955) 
and Dzhcrbashyan (1966). 

Theorem 3.1. The following results hold: 



dz 1 
d m 
dz" 



-E m (z m ) = E m (z m ) 

= 25* + £ 



fri r(l-rm/n) 



, n = 2,3, ... 



^ m— 1 

Ba(z) = — ^ Ei(z™ cxp(i27rr/m)) 



n-l 



7(1 -^,«) 



^r(l-r/n) 

w/iere 7(0, z) denotes the incomplete gamma function, defined by, 

7(0, «) - / e - ** 0-1 ^. 
Jo 



, n = 2,3, ... 



(3.1) 
(3.2) 

(3.3) 

(3.4) 

(3.5) 



In order to establish the above formulas, we observe that (3.1) and (3.2) readily follow from the definition 
(1.2). For proving the formula (3.3), we recall the identity 



m— 1 



r=0 



E r i to if k = (mod to) 

expl&fcr/m]^ (3.6) 



if k ^ (mod to). 



By virtue of the results (1.1) and (3.6), we find that 

m— 1 

E a (ze i2 ^ m ) = mE am (z m ),me N (3.7) 



r=0 

which can be written as 

m— 1 



m 

r=0 



and the result (3.3) now follows by taking a = m/n. To prove the relation (3.4), we set to = 1 in (3.1) and 
multiply it by exp(— z) to obtain 



_d 



e- z E ± (z~) =e~ z V — --. (3.9) 

y J ^ r(i - r/m) 

r— 1 v ' y 

On integrating both sides of the above equation with respect to z and using the definition of incomplete 
gamma function (3.5), we obtain the desired result (3.4). An interesting case of (3.8) is given by 

E 2a {z 2 ) = \[E a {z) + E a (-z)]. (3.10) 

4. Basic properties 

This section is based on the paper of Berberan-Santos (2005). From (1.1) and (1.2) it is not difficult to 
prove that 

E a (-x) = E 2a {x 2 ) - xE 2a . l+a {x 2 ) (4.1) 

and 

E a (-ix) = E 2x {-x 2 ) - ixE 2a . 1+a (-x 2 ), i = V^T. (4.2) 

It is shown in Berberan-Santos (2005, p. 631) that the following three equations can be used for the direct 
inversion of a function I(x) to obtain its inverse H(k): 

„c/c roc 

H(k) = — / [5ft[J(c + ioj)] cos(ku) - $s[I{c + iiv)} sm(kuj)duj (4.3) 
I" Jo 

= / SR[7(c + iuj) cos(fcw)dw, k > (4.4) 

JO 

2 e cfe roo 

= / 9f[7(c + iu)} sin(fcw)dw, k > 0. (4.5) 

7f JO 

With the help of the results (4.2) and (4.4), it yields the following formula for the inverse Laplace transform 
77(fc) of the function E a (—x). 

2 f 00 

H a (k) = - / E 2a {-t 2 ) cos(fcf)dt, k > 0,0 < a < 1. (4.6) 
t Jo 
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In particular, the following interesting results can be derived from the above result. 



2 f°° 2 f°° 

Hi(k) = - cosh(ii) cos(fci)dt = - / cos(i) cos(fct)dt = 6(k-l),i= x/^T (4.7) 
^ Jo ^ Jo 



2 /' 

£Ti (fc) = - / e~* 2 cos(/rf)di = -^e"^ (4i 
2 ^ Jo 



1 _ii 



2 r 00 2 

H?_(k) = -j e* 2 erfc(t 2 )cos(H)dt. (4.9) 

Another integral representation of H a (k) in terms of the Levy one-sided stable distribution L a (k) was given 
by Pllard (1948) in the form 

H a (k) = ^k-^ 1+ ^L a (k-i). (4.10) 
The inverse Laplace transform of E a {—x^), denoted by H@(k) with < a < 1, is obtained as 

H?(k)= r 'tfL a (t)L p (ktf)dt, (4.11) 
Jo 

where L a {t) is the one-sided Levy probability density function. From Berberan-Santos (2005, p. 432) we have 

1 f°° 

H a (k) = - [E 2a (-LJ 2 ) cos(fcw) + LoE 2a i +a {-uj 2 ) sin(fcw)]dw, < a < 1. (4.12) 

Jo 

Expanding the above equation in a power series, it gives 

1 

H a (k) = - V b n {a)k n , < a < 1 (4.13) 



7T 

n=0 



with 



C-OO 

b {a) = / £ 2a (-t 2 )dt. (4.14) 
Jo 



The Laplace transform of the equation (4.13) is the asymptotic expansion of E a {— x) as 

n=0 

5. Recurrence relations 

By virtue of the definition (1.2), the following relations are obtained in the form of 
Theorem 5.1. We have 

E a ^(z) = zE aiCt+0 (z) + (5.1) 

E aJJ (z) = pE a ^ +1 {z) + az^-E aJj+1 (z) (5.2) 

d™ 1 

- [z^E^z )] = z^- m - 1 E a>p _ m {z a ), W ~ m ) > 0, m e N. (5.3) 



dz 



— E a p(z) = . (5.4) 

az az 



The above formulae arc useful in computing the derivative of the Mittag-Leffler function E a .p(z). The 
following theorem has been established by Saxena (2002): 

Theorem 5.2. If 3?(a) > 0, Sft(/3) > and r € N then there holds the formula 

z r E a! p +ra (z)=E at p{z)-^—— -. (5.5) 

n=0 [f + > 

Proof. We have from the right side of (5.5), 

E a ,p{z) - J2 r(/3 + no) = ^ r(/3 + no) ' 

n=0 v ' n=r v ' 

Put n — r = korn = k + r. Then 



z n 00 z k+r 



^ T((3 + na) ^ T( h w + /,-n I 

n=r v ' k=0 

= Z~' E a ,f} +ra (z). 

For r = 2, 3, 4 we obtain the following corollaries: 

Corollary 5.1. If$l(a) > 0, > f/ien i/iere holds the formula 



z*E a , f3+2a (z) = E aA z) - ^ - ^^-y (5.6) 



Corollary 5.2. If$t(a) > 0, Ji(/3) > i/ien i/iere holds the formula 

I z z 2 

z 3 E aJj+3a (z) = E a ^z) ~ W) - - T{2a + f3 y ( 5 - 7 ) 

Corollary 5.3. If$l(a) > 0, 3f(/3) > £/ien iftere fto/ds £Ae formula 

1 z z^ z*^ 

z*E a , p+4a {z) = E aJj (z) - — - j^-^ - r(2a + /3) - r( 3a + /3)- (5 ' 8) 

Remark 5.1. For a generalization of the result (5.5), see Saxena, Kalla and Kiryakova (2003). 

6. Asymptotic expansions 

The asymptotic behavior of Mittag-LcfHcr functions plays a very important role in the interpretation of 
the solution of various problems of physics connected with fractional reaction, fractional relaxation, fractional 
diffusion and fractional reaction-diffusion etc in complex systems. The asymptotic expansion of E a {z) is 
based on the integral representation of the Mittag-Leffler function in the form 

Ea{z) = idb ta ~t<* -*f } dt ' jR(a) > °' a,z e c ' (6J) 

where the path of integration f2 is a loop starting and ending at — oo and encircling the circular disk \t\ < \z\i 
in the positive sense, | argt| < it on O. The integrand has a branch point at t = 0. The complex i-plane is 
cut along the negative real axis and in the cut plane the integrand is single- valued, the principal branch of 
t a is taken in the cut plane. (6.1) can be proved by expanding the integrand in powers of t and integrating 
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term by term by making use of the well-known Hankcl's integral for the reciprocal of the gamma function, 
namely 



1 



T(f3) 2m J Ha 



■X 



rdC- 



(6.2) 



The integral representation (6.1) can be used to obtain the asymptotic expansion of the Mittag-LefHer 
function at infinity (Erdclyi, et at., 1955). Accordingly, the following cases are mentioned below: (i): If 
< a < 2 and \i is a real number such that 



na • r i 

— < /i < min[7T, na\, 
then for TV* £ N, N* ^ 1 there holds the following asymptotic expansion: 

N* 



E a {z) = ±zW°ev>{zi)-'E m 
as \z\ — > oo, | argz| < fi; and 



1 1 



J r(l — ar) z r 



O 



1 



N* 



1 1 



r—l v ' 



o 



y JV* + l 



as \z\ — > oo, \l < | argz| < ir. (ii) : When a > 2 then there holds the following asymptotic expansion: 



N* 

_ , . 1 V"^ i r / 2t17TZ . i ^—v 1 1 

#«(z) = - > £™ expexp )Z"\ - > — r— + 

a ^— ' a ^ r(l — ar) z r 



O 



1 



as \z\ — > oo, | argz| < and where the first sum is taken over all integers n such that 



arg(z) + 2nn\ < 



air 



(6.3) 



(6.4) 



(6.5) 



(6.6) 



(6.7) 



The asymptotic expansion of E a ^(z) is based on the integral representation of the Mittag-Lefner function 
E a ,[){z) in the form 



E 



a,/3 



1 f t a -Pexp(t) 



dt, Sft(a) > 0, 9?(/3) > 0, z, a, (3 £ C, 



(6.8) 



la t a - z 

which is an extension of (6.1) with the same path. As in the previous case, the Mittag-LefHer function has 
the following asymptotic estimates: (Hi): If < a < 2 and /x is a real number such that 



ira 

— < fj, < min[7r, na\ 



then there holds the following asymptotic expansion: 



AT* 



EaA*) = ~z (1 - p)/a exp(z* ) - Y, WT ^ + ° 
a J T(p — ar) z r 



1 



z N» + l 



(6.9) 



(6.10) 



as \z\ — > oo, | argz| < /x; and 



N* 



E a ,p{z) = -Y i - 



1 1 



+ 



1 



y N' + l 



^ T(0 - ar) z r 

r—l v 7 

as \z\ — > oo, /j, < I argz| < n. (iv): When a > 2 then there holds the following asymptotic expansion: 

JV* 

> z« exp|cxp( )z°| ^ — 



= exp[exp(^)^]^ - £ 



1 1 



^T(/3-ar) zr 



O 



1 



(6-11) 



(6.12) 



as \z\ — > oo, | argz| < ^ and where the first sum is taken over all integers n such that 



arg(z) + 2irn\ < 



(6.13) 



7. Integral representations 

In this section several integrals associated with Mittag-LefRcr functions are presented, which can be 
easily established by the application by means of beta and gamma function formulas and other techniques, 
see Erdelyi, et al. (1955), Gorenflo ct al. (1997, 2002). 



/ c-<E a (Cz)d( = - , \z\ <l,a€C, K(a) > 

f°° 1 

/ e-'xP^Ea {x a z)dx = , \z\ < 1, a, G C, R(a) > 0, R(/3) > 

Jo 1 - z 

(x - O^^OdC - T{P)x fi E a , p+1 {x a ), > 



t,a — l 



f°° i \ rnlq"-/ 3 



(s Q T a) 



where a, /? € C and 



(7.1) 
(7.2) 
(7.3) 
(7.4) 
(7.5) 



dz' 
2 



;E a A z ) 



E a (—x a ) = — sin(a7r 



l2) L T 



C"- 1 cos(xC) 



A „-dC, a G C, 3?(a) > 
./n l + 2C Q cos(o:7r/2)+C 2Q 



^ r 00 /-a— 1 

E a (—x) = — sin(a7r) / — ; — , 

' ti Jo 1 + 2( a cos(aTr) + C 2q 

i x« r 00 

E a (—x) = 1 1 / arctan 

2a 7T J 



e _Cx °C dC« G C,R(a) > 



C" + cos(a7r) 
sin(a7r) 



e" Ca:Q C d(,a G C,R(a) > 0. 



(7.6) 
(7.7) 
(7.8) 



Note 7.1. Equation (7.7) can be employed to compute the numerical coefficients of the leading term of the 
asymptotic expansion of E a (—x). Equation (7.8) yields 



a f°° 
b (a) = -r(a) sin(2a7r) / - 
n Jo 1 



+ 2C 2Q cos(2a7r) + C 4adC ' a < 2' 



From Berberan-Santos (2005) and Gorenflo et al.(1997) the following results hold: 



(7.9) 



E a {-x) = — r ^-p-dt,0<a<l,aeC. 

TT J X 2 + t 1 



In particular, the following cases are of importance 



„ . . 2x f°° coshtit) , , , 



-*>-?/ 



TT Jo a; 2 + t 
2a; p e* 2 erfc(-t 2 ) 
x 2 + £ 2 



dt. 



(7.10) 

(7.11) 
(7.12) 
(7.13) 
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Note 7.2. Some new properties of the Mittag-Lemer functions are recently obtained by Gupta and Dcbnath 
(2007). 



8. The H-function and its special cases 

The H-function is defined by means of a Mcllin-Barnes type integral in the following manner (Mathai 
and Saxena, 1978): 



H p,q \ Z ) - H p 



Q 





(dp , Ap ) 






(oi.Ai),. 


.,(a p ,A p 


z 




- B p , q 


z 












(6i,Bi),- 





where i — (— 1) and 



Q(s) 



{UT=i rfe + sBj)} {UU r (! - a _j g ^)} 
{UUm+i r(i - h sBj)} {u p 1=n+ i r (% + sA j)} ' 



(8.1) 



1.2) 



and an empty product is interpreted as unity; m, n,p,q € Ao with < n < p,l < m < q, Ai, Bj e a^, 6j e 
C,i = 1, ...,p;j = 1, g such that 



Ai(bj + k) ^ Bj(ai - A - 1), fc, A e N ; i=l, ...,n; j = 1, ...,m, 



(8.3) 



where we employ the usual notations: A^o = 0, 1, 2, R = (— oo, oo), i?+ = (0, oo), and C being the complex 
number field. The contour O is the infinite contour which separates all the poles of T(bj + sBj),j = 1, m 
from all the poles of T(l — ai + sAi),i = 1, ...,n. The contour could be ft = L-oe or ft = L +OQ or £1 = Li 100 
where L_oo is a loop starting at — oo encircling all the poles of T(bj + sBj),j = 1, m and ending at — oo. 
L +oc is a loop starting at +oo, encircling all the poles of T(l — — sAi), i = 1, n and ending at +oo. Li 100 
is the infinite semicircle starting at 7 — zoo and going to 7 + zoo. A detailed and comprehensive account of 
the H-function is available from the monographs of Mathai and Saxena (1978), Prudnikov et al. (1990) and 
Kilbas and Saigo (2004). The relation connecting the Wright's function p ip q (z) and the H-function is given 
for the first time in the monograph of Mathai and Saxena (1978, p. 11, eq(1.7.8)) as 



O^q 



(ai, Ax), (a p ,A p ) , 
(6i,Bi),...,(6„B,) |x 



(l-oi.Ai),...,(l-a p ,A p ) 
(0,l),(l-6i,Bi),...,(l-6 5 ,B 5 ) 



(8.4) 



where p tp q (z) is the Wright's generalized hypergeometric function (Wright, 1935, 1940); also see Erdelyi, et 
al. (1953, Section 4.1), defined by means of the series representation in the form 



p%( z ) = E 



{lX P =ir(a, + V)} z r 



where z € C, a i7 bj e C, Ai, Bj e R+, Ai =/= 0, Bj ^ 0;i = 1, ...,p;j = 1, q, 



1.5) 



9 P 



E^-E^ = A> - L 

J=l »=1 

The Mcllin-Barnes contour integral for the generalized Wright function is given by 



plpq 



(a,p, Ap) 
(b q ,B q ) 1 



_1_ f 

2™ J n UU 



=1 rfe - 



(-z)- s ds 



(8.6) 



where the path of integration separates all the poles of T(s) at the points s = — v, v £ No lying to the left 
and all the poles of ^( a j ~~ s A;')>i = 1' -■■■>P at the points s = (Aj + Vj)jAj,Vj e A^o, j = 1, ...,p lying 
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to the right. If O = (7 — 100,7 + *°°) then the above representation is valid if cither of the conditions are 
satisfied: 



(0 



A<l,|arg(-z)|< (1 o A) V ^0, 



(ii) 



A = 1, (1 + A) 7 + \ < H(<f), arg(-z) = 0, z ? 0, 5 = fl b > ~ £ a o + 



3 = 1 



3 = 1 



This result was proved by Kilbas, Saigo and Trujillo (2002). 



The generalized Wright function includes many special functions besides the Mittag-LefHer functions 
defined by the equations (1.1) and (1.2). It is interesting to observe that for A, = Bj = l,i = 1, ...,p;j — 
1, ...,q, (8.5) reduces to a generalized hypergeometric function p F q {z). Thus 



plpq 



K, 1) 



n?=ir(6;) 



pFq( a i, ■••) a p] bi, ...,b g ; z) 



(8.7) 



where aj ^ —v,j = 1, ...,p, v = 0, 1, &j 7^ — A, j = 1, q, A = 0, 1, p < q or p = q + 1, \z\ < 1. Wright 
(1933) introduced a special case of (8.5) in the form 



(a, 6; z) = 0^1 



(6, a 



I* 



1 



^ T(ar + 6) r! 

r— 



(8.8) 



which widely occurs in problems of fractional diffusion. It has been shown by Saxena, Mathai and Haubold 
(2004a), also see Kiryakova (1994), that 



(1.1) 



(M 



E a ,p(z) = iipi 
If we further take (i = 1 in (8.9) we find that 

E aA (z) = E a (z) = iV>i 

where a e C, SR(a) > 0. 



I* 



#1,2 



(0,1) 

(0,l),(l-/3,a). 



(1,1) 

(l,a) 



I* 



if 



1.2 



(0,1) 

(0,1), (0,a) 



(8.9) 



(8.10) 



Remark 8.1. A series of papers are devoted to the application of the Wright function in partial differential 
equation of fractional order extending the classical diffusion and wave equations. Mainardi (1997) has 
obtained the result for a fractional diffusion wave equation in terms of the fractional Green function involving 
the Wright function. The scale- variant solutions of some partial differential equations of fractional order were 
obtained in terms of special cases of the generalized Wright function by Buckwar and Luchko (1998) and 
Luchko and Gorenfio (1998). 

9. Mellin-Barnes integrals for Mittag-LefHer functions 

These integrals can be obtained from the identities (8.9) and (8.10). 
Lemma 9.1. If$l(a) > 0,3?(/3) > and z e C the following representations are obtained: 



1 T(s)T(l-s) 



E a .A z ) 



r(i (\s 
j_ f t+tca r(a)r(i - s) 

2^i 



I 



7 - JO o r(/3 - as) 
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(-z)- s ds 



(9.1) 
(9.2) 



where the path of integration separates all the poles of T(s) at the points s = —is, v = 0, 1, ... from those of 
r(l — s) at the points s = 1 + v, v = 0, 1, .... 

On evaluating the residues at the poles of the gamma function T(l — s) we obtain the following analytic 
continuation formulas for the Mittag-Lefner functions: 

Ea{z) ~ ™ 7 7 - ioo m- a s) { - z) ds = -^f(i-^ (9 - 3) 



and 



i n +i °° r( s )r(i - 8 ) 



-k 



-55JU ^T ( -= r ' ds -£n^SE5- (9 - 4) 

10. Relation with Riemann-Liouville fractional calculus operators 

In this section, we present the relations of Mittag-Lefner functions with the left and right-sided operators 
of Riemann-Liouville fractional calculus, which are defined below. 



(Jo + 0)O»O = pn *) Q_ V(t)dt,K(a) > (10.1) 
r (") Jo 

1 Z" 00 

= fR 7 (< " s)"" V(*)dt, »(«) > ( 10 - 2 ) 

(W)(z) = (-^) [Q %-~ {Q V](Z) 

= f(r3M)(-d^J I (*-*)- W ^)d*,»(a)>0 (10.4) 

where [a] means the maximal integer not exceeding a and {a} is the fractional part of a. 

Note 10.1. The fractional integrals (10.1) and (10.2) are connected by the relation (Kilbas, 2005, p. 118) 

[i«4>(i/t)}( X ) = x a - 1 (5 + r Q -vw]) (I) . 



,a^0 (10.5) 



Theorem 10.1. Let 3?(a) > and 3?(/3) > then there holds the formulas 

r /3-l 



{I^[t^E a .,(an\)(x) ' 



E a .j3(ax a ) 



a 

(IZ + [E a (at a )])(x) = - a [E a (ax a ) - 1) ,a ? (10.6) 

w/wc/i 6y virtue of the definitions (1.1) and (1.2) can be written as 

(I^y-'E^iaeMx) = x a+ P- l E a , a+ p(ax a ) (10.7) 
(I£ + [E a (at a )})(x) = x a E a , a+1 (ax a ). (10.8) 
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Theorem 10.2. Let 5t(a) > and 3?(/3) > then there holds the formulas 



(I«[t- a -PE a!0 (at- a )])(x) 



E a ^(ax a ) 



(I*[ t -<x-i-E a (at- a )})(x) = -x a - 1 [E a (ax- a ) - l],a ^ 0. 



(10.9) 
(10.10) 



Theorem 10.3. Let < $l(a) < 1 and > $l(a) then there holds the formulas 

Ji-a-\ 



{D^[t^E aS {at a )\)(x) 



T(P-a) 



+ ax f3 - 1 E a .j 3 (ax a ) 



(Dg + [E a (at a )})(x) = ^— - + aE a (ax a ). 



(10.11) 
(10.12) 



Theorem 10.4. Let $l(a) > and 3?(/3) > [9?(a)] + 1 then there holds the formula 



(Dl[t a -^E a ^(at- a )])(x) = 



-0 



T(J3-a) 

11. Generalized Mittag-Leffler type functions 



ax- a ~ fi E a .p{ax- a ). 



(10.13) 



By means of the series representation a generalization of (1.1) and (1.2) is introduced by Prabhakar 
(1971) as 



(7), 



n=0 



\T(/3 + an) 



,a,p,je C,$(a) > 0,3?(/3) > 0, 



(11.1) 



where 



(7)„=7(7+l)-(7 + n-l) = ^±^ 



whenever T(-f) is defined, (7)0 = 1,7 ^ 0. It is an entire function of order p — [3?(a)] _1 and type a = 
- [{^(a)} 3 ^"'] P . It is a special case of Wright's generalized hypergeometric function, Wright (1934,1935) 
as well as the H- function (Mathai and Saxena, 1978). For various properties of this function with applications, 
see Prabhakar (1971). Some special cases of this function are enumerated below. 



E a (z) = Ei A (z). 
E a , p (z) = El (z). 

*iK + f hz) = (i + «7 - P)Kd z ) + E le-^)- 

cj>(a,f3;z)=r(f3)E^(z) 



(11.2) 
(11.3) 
(11.4) 
(11.5) 



where <j){ot, [3; z) is the Rummer's confluent hypergeometric function. Ej, {z) has the following representa- 
tions in terms of the Wright's function and H-function. 



(7,1). 

(/3, a)" 



1,1 





(1-74) 








(0,1), (1-/3, a)_ 



27rwr( 7 ) 



r( s )r( 7 - s ) 

T((3-as) 
12 



(-zy s ds,^(j) > 



(11.6) 
(11.7) 
(11.8) 



where iV'i(') and H^i') are respectively Wright generalized hypergeometric function and the H-function. In 
the Mellin-Barnes integral representation, uj = V^T and the c in the contour is such that < c < 3?(7) and 
it is assumed that the poles of T(s) and T(7 — s) are separated by the contour. The following two theorems 
are given by Kilbas, Saigo and Saxena (2004). 

Theorem 11.1. If a, (3, 7, a £ C, 3?(a) > 0, 3?(/3) > n, > then for n e N the following results hold: 



In particular, 



and 



J 71 



J-n-lp7 /„ 7 t>\ 
^a,fi-n\ aZ I 



dz 



d 

dz 



^[zf-^ftaz)] 



Ji-n-l 



<t>(r,P- n;az). 



(11.9) 
(11.10) 
(11.11) 



Theorem 11.2. If a, (3, 7, a, v, a £ C, Sft(a) > 0, &(/?) > 0, ^(7) > 0, SR(i/) > 0, SR(ct) > f/ien 
"(a: - ^"^X* ~ t) a ]f- 1 ^,,(a* a )dt = x^E^^x"). 



J 

Jo 



The proof of (11.12) can be developed with the help of the Laplace transform formula 



P- l El t0 {ax a ) (s)=s-P(l-as- a ) 



a\-7 



(11.12) 



(11.13) 



where a,/?, 7, a € C, Sft(a) > 0,5R(/3) > 0,3?( 7 ) > 0, s e C,3?(s) > 0,|as~ Q | < 1. For 7 = 1, (11.13) reduces 
to 



L [x f ^ 1 E a .p{ax a )} (s) = s" /3 (l - as-")- 1 . 
Generalization of the above two results is given by Saxena (2002). 



L [I?- 1 £%„(<*<>)] (s) 



T(6) 



(6,1), (p, a) a_ 

(7,/3) 's Q 



(11.14) 



(11.15) 



where Sft(/3) > 0, 3?( 7 ) > 0, 5R(s) > 0, &(p) > 0, s > \a\ ^ , $1(6) > 0. 

Relations connecting the function defined by (11.1) and the Riemann-Liouvillc fractional integrals and 
derivatives are given by Saxena and Saigo (2005) in the form of nine theorems. Some of the interesting 
theorems are given below. 

Theorem 11.3. Let a>0, /3>0,7>0 and a € R. Let Iff be the left-sided operator of Riemann-Liouville 
fractional integral. Then there holds the formula 



(Iff^El^Mx) = X ^-'El a+1 (axP). 



(11.16) 



Theorem 11.4. Let a>0,/3>0,7>0 and a G R. Let I™ be the right-sided operator of Riemann-Liouville 
fractional integral. Then there holds the formula 



(I«[t- a ^E s 0tJ (at-P)})(x) = x-itf^iax- 



(11.17) 
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Theorem 11.5. Let a>0,/3>0,7>0 and a £ R. Let D$ + be the left-sided operator of Riemann-Liouville 
fractional derivative. Then there holds the formula 

{D^-'E$^))){x) = x'— l Ef^_ a (ax^. (11.18) 



Theorem 11.6. Let a > 0,(3 > 0,7 — a + {a} > 1 and a £ R. Let be the right-sided operator of 
Riemann-Liouville fractional derivative. Then there holds the formula 

(D a _[t^E s ^(ar^})(x) = x^E s ^_ a (ax-^. (11.19) 

In a series of papers by Luchko and Yakubovich (1990, 1994), Luckho and Srivastava (1995), Al-Bassam and 
Luchko (1995), Hadid and Luchko (1996), Gorcnflo and Luchko (1997), Gorenflo et al. (2000), Luchko and 
Gorenflo (1999), the operational method was developed to solve in closed forms certain classes of differential 
equations of fractional order and also integral equations. Solutions of the equations and problems considered 
are obtained in terms of generalized Mittag-Leffier functions. The exact solution of certain differential 
equation of fractional order is given by Luchko and Srivastava (1995) in terms of the function (11.1) by using 
operational method. In other papers, the solutions are established in terms of the following functions of 
Mittag-Leffier type: If z,p,f3j £ C,$l(aj) > 0, j = 1, ...,m and m £ N then 

JW«.i,A)W»» = Sn^E+ftJS- (1L20) 

For m = 1, (11.20) reduces to (11.1). The Mcllin-Barnes integral for this function is given by 



1 p +lo ° T(s)T(p-a) 

27TlL(p) J 1 _ loo r{(3j - CLjS) 



(—z) s ds,i 



r(p) 



(i-p,i) 

(0,1), (1-Pj, aj ),j = l,...,m 



(11.21) 
(11.22) 



where 0,7 < 3?(/o),9?(p) > and the contour separates the poles of T(s) from those of T(p — s). $t(aj) > 
Q,j = 1, m, arg(— z) < ir. The Laplace transform of the function defined by (11.20) is given by 



L[E p (( aj , /3j) 1>m ; -t)](s) 



sT(p) 



(P,l),(l,l) 

{aj,Pj)i, m 



(11.23) 



where SR(s) > 0. 



Remark 11.1. In a recent paper, Kilbas, Saigo and Saxena (2002) obtained a closed form solution of a 
fractional generalization of a free electron equation of the form: 



D>(r) = A f t s a(r - t)E b p 4+1 (ii/f)dt + P^Ej a+1 (ii/f ), < r < 1, ; 
Jo 



'-I 



(11.24) 



where 6, A <G C,u,f3 £ R +1 a > 0, p > 0;a > —l,p > —1,5 > —1 and E h p 5+1 {-) is the generalized Mittag- 
Leffier function given by (11.1), and a(r) is the unknown function to be determined. 

Remark 11.2. The solution of fractional differential equations by the operational methods are also obtained 
in terms of certain multivariate Mittag-Leffier functions defined below: The multivariate Mittag-Leffier 
function of n complex variables z\,...,z n with complex parameters a\, a n , b £ C is defined as 



oo L\-\r... J rL n —k 

E (ai ^ an)ib (z 1 ,...,z n ) = 

k=0 L!,...,L„>0 



L u ...,L n J T(b±Y^ =1 ajLj) 



(11.25) 
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in terms of the multinomial coefficients 



k \ fc' 

t t = t \ 't , ,k,Lj,eN ,j = l,...,m. (11.26) 

Another generalization of the Mittag-LefHer function (1.2) was introduced by Kilbas and Saigo (1995) 
in terms of a special function of the form 



fc-i 



,,,,(*) = co = 1, c k = n J}?}!™ 1^ . ! < £ v " = {0.:i 2....} (1-1.27) 



fc=0 



i=0 



T(a(im + /3 + 1) + 1) : 



where an empty product is to be interpreted as unity; a, e C are complex numbers and m G R 1 di(a) > 
0, m > 0, ct(im + 0) £ Z~ = {0, — 1, —2, ...}, i = 0, 1, 2, ... and for m = 1 the above defined function reduces 
to a constant multiple of the Mittag-Lcfner function, namely 



E a ,i,p{z) - T{a0 + l)E a . aP+1 (z), 



(11.28) 



where 5R(a) > and a(i + 0) ^ It is an entire function of z of order [3?(a)] _1 and type a = 1/m, 

see Gorenflo et al. (1998). Certain properties of this function associated with Riemann-Liouville fractional 
integrals and derivatives are obtained and exact solutions of certain integral equations of Abel-Volterra type 
are derived by their applications (Kilbas and Saigo, 1995, 1996). Its recurrence relations, connection with 
hypergeometric functions and differential formulas are obtained by Gorenflo, Kilbas and Rogosin (1998), also 
see, Gorenflo and Mainardi (1996). In order to present the applications of Mittag-LefHer functions we give 
definitions of Laplace transform, Fourier transform, Riemann-Liouville fractional calculus operators, Caputo 
operator and Weyl fractional operators in the next section. 

12. Laplace and Fourier transforms, fractional calculus operators 

We will need the definitions of the well-known Laplace and Fourier transforms of a function N(x, t) and 
their inverses, which are useful in deriving the solution of fractional differential equations governing certain 
physical problems. The Laplace transform of a function N(x,t) with respect to t is defined by 



/>oo 

L{N(x,t)} = N(x,s) = / c- st N(x,t)dt,t > 0,x £ R 
Jo 

where ^St(s) > and its inverse transform with respect to s is given by 

L-^Nix, a)} = L~ 1 {N(x, s); t} = — / e st N(x, s)ds. 
The Fourier transform of a function N(x, t) with respect to x is defined by 

/oo 
e lkx N(x, t)dx. 

The inverse Fourier transform with respect to k is given by the formula 



1 f°° 

F-^F^M)} = 7^ / c' lkx F*{k,t)dk. 

271" J-oo 



(12.1) 



(12.2) 



(12.3) 



(12.4) 



From Mathai and Saxena (1978) and Prudnikov et al.(1990, p. 355, eq(2.25.3)) it follows that the Laplace 
transform of the H-function is given by 







(dp )-<4p) 




L{tP^H^ q n 


zt a 




} = s 









— b n p+l,q 





(1— p,cr) ,(a p ,A P )' 


zs- a 






(b q ,B g ) 



(12.5) 
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where Sft(s) > 0, + crmini< J < m (-^-) > 0, cr > 0, 

p m q 

(12.6) 

By virtue of the cancelation law for the H-function (Mathai and Saxena, 1978) it can be readily seen that 

(12.7) 



^ n p m q 

| argz| < -nCt, Q, > 0,0 = A, - ^ A i + ^2 B j- ^ B r 

z— 1 i= n+1 j=l j=m+l 













(a p ,A p ),(p,a)~ 


L-^s-'HZf 


zs a 


(b,>S,)_ 


J n p+l,q 




{b q ,B q ) 



where a > 0, 3?(s) > 0, + (rmaxi< i < n (^ i )] > 0, | argz| < ^irQi, Oi > 0, Oi = O — p. In view of the 
results 



(12.8) 



J_ i (x) = \ COS X 

2 V 7TX 

the cosine transform of the H-function (Mathai and Saxena, 1978, p. 49) is given by 



p CO 



a 



(p,,,),(a-a p ,A p ).X^^), 



(12.9) 



where + ^mini< J -< m (^ r )] > 0, 3R[p + maxi<j< n (^-i) < 0, |arga| < |7rQ,Q > 0; O = Y!j=i B j ~ 
EUn+i B J + EU A i n=„+i ^> and fc > 0. 

The definitions of fractional integrals used in the analysis arc defined below. The Riemann-Liouville 
fractional integral of order v is defined by (Miller and Ross, 1993, p. 45) 



{t-u) v ~ 1 f(x,u)du, 



(12.10) 



where 5i(z/) > 0. Following Samko, Kilbas and Marichev (1993, p. 37) we define the Riemann-Liouville 
fractional derivative for a > in the form 



oD?f(x,t) 



1 



(t-u) n - a -\f(x,u)du, n=[a] + l 



r(n - a) dt n J 

where [a] means the integral part of the number a. From Erdelyi, et al. (1954b, p. 182) we have 



(12.11) 



L{ D^f{x,t)} = s- v F{x,s), (12.12) 

where F(x, s) is the Laplace transform of f(x, t) with respect to t, 3?(s) > 0, 3?(f) > 0. The Laplace transform 
of the fractional derivative defined by (12.11) is given by Oldham and Spanier (1974, p.l34,eq (8.1.3)) 

n 

L{ Q D<?f{x,t)} = s a F{x,s)-Y,s k ~ 1 zDr k f{x,t)\ t=0 , n-l<a<n, (12.13) 

fc=i 

In certain boundary-value problems arising in the theory of visco-elasticity and in the hereditary solid 
mechanics the following fractional derivative of order a > is introduced by Caputo (1969) in the form 



D?f(x,t) 



, (t-T) m - a - 1 f (m \x,t)dt,m-l<a<m,n(a) >0,me N (12.14) 
T(m- a) J Q 

d m 

f(x, t), if a = m, 



dt 



(12.15) 



where Jpr/ is the m-th partial derivative of the function f(x,t) with respect to t. The Laplace transform 
of this derivative is given by Podlubny (1999) in the form 
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m—1 



L{D?f(t); s} = s a F(s) - J2 ■s a - k - 1 .f ik) (0+), m - 1< a < m. (12.16) 



k=0 



The above formula is very useful in deriving the solution of diffcrintegral equations of fractional order 
governing certain physical problems of reaction and diffusion. Making use of the definitions (12.10) and 
(12.11) it readily follows that for f(t) = t p we obtain 

oort p = r^+^+i) *^'^ > °' ^ > -i; * > o (12.17) 

and 

DW = J {p+1 ^, tr-»W) < 0, K(P) > -1; t > 0. (12.18) 
On taking p = in (12.18) we find that 

D» t [l]= ^_ t-",t>0&{y)<l. (12.19) 

From the above result, we infer that the Ricmann-Liouville derivative of unity is not zero. We also need the 
Weyl fractional operator defined by 

1 d™ /"* 

where n = [p] + 1 is the integer part of p > 0. Its Fourier transform is (Metzler and Klaftcr, 2000, p. 59, 
A.ll) 

F{^D^f(x)} = (ikrf*(k), (12.21) 
where we define the Fourier transform as 

/oo 
h(x)exp(iqx)dx. (12.22) 
-oo 

Following the convention initiated by Compte (1996) we suppress the imaginary unit in Fourier space by 
adopting a slightly modified form of the above result in our investigations (Metzler and Klafter, 2000, p. 59, 
A.12). 

F{- 00 D£f(x)} = -\k\'>r(k), (12.23) 

instead of (12.21). 

We now proceed to discuss the various applications of Mittag-Leffler functions in applied sciences. In 
order to discuss the application of Mittag-Leffler function in kinetic equations, we derive the solution of two 
kinetic equations in the next section. 

13. Application in kinetic equations 

Theorem 13.1. J/3?(i/) > then the solution of the integral equation 

N(t) -N = -c\Di v N(t) (13.1) 

is given by 

N{t)=N E u {-c v t v ) (13.2) 
where E v {t) is the Mittag-Leffler function defined in (1.1). 
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Proof. Applying Laplace transform to both sides of (13.1) and using (12.12) 

N(s) = L{N(t);s} = Nos-^l + (a/c) - "] -1 - 

By virtue of the relation 



it 



gives 



(13.3) 



fc=0 

j-i/k 



N J2(- 1 ) kcVk r{1 + l/k) = N oEA-^n- (13-6) 



L-H*-"} = fyj&iP) > °> ^( s ) > 0, s e C, (13.4) 

it is seen that 

oo 

L-^Nos-^l + (s/c)-"}- 1 } = N Y,( k -l) h <? k L- 1 {8-' k - 1 } (13.5) 

t v 

- irr "f? 

k=0 

This completes the proof of Theorem 13.1. 

Remark 13.1. If we apply the operator D" from the left to (13.1) and make use of the formula 

oD v t [\]= 1 f~V >(),%(»/) <1 (13.7) 
we obtain the fractional diffusion equation 

D»N(t) - N v T ^ V _ v) = -c u N{t),t > 0, < 1 (13.8) 
whose solution is also given by (13.6). 

Remark 13.2. We note that Haubold and Mathai (2000) have given the solution of (13.1) in terms of the 
series given by (13.5). The solution in terms of the Mittag-Leffler function is given in Saxena, Mathai and 
Haubold (2002). 

Alternate procedure. We now present an alternate method similar to that followed by Al-Saqabi and 
Tuan (2006) for solving some diffcrintegral equations, also, see Saxena and Kalla (2008) for details. 

Applying the operator {-c v ) m D^ mv to both sides of (13.1) we find that 

{-c v ) m n D~ mv N{t) - {-c v ) m+1 D- v(m+1) N{t) = N Q D^ m ''[i\,m = 0, 1,2, ... (13.9) 
Summing up the above expression with respect to m from to oo, it gives 

oo oo oo 

]T (-cToDt m »N(t) - (-c v ) m+1 D; v{m+1) N(t) = NoY, {-c v ) m *Di mv [l\ 

m— m— m— 

which can be written as 

oo oo oo 

m—0 m—1 m—0 

Simplifying the above equation by using the result 



D7 v t"- 1 = r ^ j p+v-i (13.10) 

r(M + ^) 

where min{3?(^), > 0, we obtain 
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00 j-mv 

N(t) = N J2 { -cn- WT ^- (13.11) 

m— 

for m = 0,1,2,... Rewriting the series on the right in terms of the Mittag-Leffler function, it yields the 
desired result (13.6). The next theorem can be proved in a similar manner. 

Theorem 13.2. If min{3?(z/), Jt(/x)} > then the solution of the integral equation 

N(t) - iVot^ 1 = -c\Dt V N{t) (13.12) 

is given by 

N(t) = N T(n)f i - 1 E v>li (-^), (13.13) 

where E u ^(t) is the generalized Mittag-Leffler function defined in (1.2). 

Proof. Applying Laplace transform to both sides of (13.12) and using (1.11), it gives 

N{s) = {N(t);s} = 7V r( M ) S ^[l + {s/c)-"]- 1 . (13.14) 

Using the relation (13.4), it is seen that 

oo 

L-^NoT^s-^l + (s/c)-}- 1 } = JV £(-l) fc c wfc L- 1 {s-"-' fc } 



k=0 

°° f vk 



k) 



NoT^-'E^i-c^). (13.15) 



This completes the proof of Theorem 13.2. 



Alternate procedure. We now give an alternate method similar to that followed by Al-Saqabi and Tuan 
(2006) for solving the diffcrintegral equations. Applying the operator {—c' J ) m D^ mlJ to both sides of (13.12), 
we find that 

{-c v ) m Dt mv N{t) - {-c v ) m+1 D- v{m+1) N{t) = Noi-cToDt™^- 1 (13.16) 
for m = 0, 1, 2, .... Summing up the above expression with respect to m from to oo, it gives 

oo oo 

J2 {-c u ) m Q D- mu N{t) - (-c v ) m+1 D- v{m+1) N(t) 

ni—0 m— 

oo 

= N J2 {-d') m Dt mv lP- x 

which can be written as 

oo oo oo 

(-c u ) m Q D- mv N(t) - (-c v ) m Dj mv N(t) = N Y {-c u ) m Q D- mu t^-\ 

m—0 m— 1 m=0 

Simplifying by using the result (13.10) we obtain 

N(t) = N T(p) V {-c»T- — -; m = 0, 1, 2, ... (13.17) 

19 



Rewriting the series on the right of (13.17) in terms of the generalized Mittag-LefHer function, it yields the 
desired result (13.5). Next we present a general theorem given by Saxena et al.(2004). 



Theorem 13.3. If c > 0, 3?(^) > then for the solution of the integral equation 

N(t) - N f(t) = -c%D- v N{t) 
where f(t) is any integrable function on the finite interval [0,6], there exists the formula 

(-i/M) 



N(t) 



= cN [ H\ l 2 
Jo 



c»(t- T y 



/(r)dr, 



(13.18) 



(13.19) 



where il 1 '^-) is the H-function defined by (8.1). 



The proof can be developed by identifying the Laplace transform of N () + c v $D^ V N (t) as an H-function 
and then using the convolution property for the Laplace transform. In what follows, Ep (■) will be employed 
to denote the generalized Mittag-LefHer function, defined by (11.1). 

Note 13.1. For an alternate derivation of this theorem see Saxena and Kalla (2008). 
Next we will discuss time-fractional diffusion. 

14. Application to time-fractional diffusion 

Theorem 14.1. Consider the following time-fractional diffusion equation 



^N{x, t) = D^N(x, t), < a < 1, 



(14.1) 



where D is the diffusion constant and N(x, t = 0) = 8(x) is the Dirac delta function and lim x ^± 00 N(x, t) 
0. Then its fundamental solution is given by 



N{x,t) = ±rHl$ 





(!,«)" 


[M 2 


Dt a 


(1,2). 



(14.2) 



Proof. In order to find a closed form representation of the solution in terms of the H-function, we use the 
method of joint Laplace- Fourier transform, defined by 



/>oo />OC 



N*(k,s) = 



-st-\-ikx 



N(x, t)dx dt 



(14.3) 



where, according to the convention followed, ~ will denote the Laplace transform and * the Fourier transform. 
Applying the Laplace transform with respect to time variable t, Fourier transform with respect to space 
variable x and using the given condition, we find that 



s a N*(k, s) - s"- 1 = -Dk 2 N*(k, s) 

and then 



N*(k lS ) 



s a-l 



s a + Dk 2 ' 

Inverting the Laplace transform, it yields 

N * {k > t] = £ " 1{ ^TW } = E ^- DkHa )> ( 14 - 4 ) 

where E a {-) is the Mittag-LefHer function defined by (1.1). In order to invert the Fourier transform, we will 
make use of the integral 
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f 

Jo 



cos(kt)E a ^(-at 2 )dt = 



'k 2 


(/?>«)" 


a 


(1,2). 



where SK(a) > 0, &(/?) > 0, k > 0, a > 0; and the formula 



1 Z" 00 1 
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e~ ikx f(k)dk=- I f(k)cos(kx)dk. 



Then it yields the required solution as 



1 



M 2 ,(i,i), a,*), (i,i) 



£>£a 1(1,2), (1,1), (1,1) 

Note 14.1. For a = 1, (14.7) reduces to the Gaussian density 

JV(x,t) 



-if 



i.i 



|x| 2 ,(l, a )' 



£>*al(1.2) 



(14.5) 



(14.6) 



(14.7) 



1 \x\ z 



(14.8 



Fractional space-diffusion will be discussed in the next section. 

15. Application to fractional-space diffusion 

Theorem 15.1. Consider the following fractional space- diffusion equation 

§- t N(x,t) = D-^N(x,t),0<a<l, (15.1) 

where D is the diffusion constant, is the operator defined by (12.20) and N(x, t = 0) = S(x) is the Dirac 
delta function and Ym\ x ^± 00 N(x, t) = 0. Then its fundamental solution is given by 

i r i-i (i,^)>(i-3)i 

N(x,t) = — H% -LL . (15.2) 

a\x\ ^ [{Dt)~ (i,i), (i,i)_ 

The proof can be developed on similar lines to that of the theorem of the preceding section. 

16. Application to fractional reaction-diffusion model 

In the same way, we can establish the following theorem, which gives the fundamental solution of the 
reaction-diffusion model given below. 

Theorem 16.1. Consider the following reaction- diffusion model 



N(x,t)=r] -coD*N{x,t),0<f3<l 



(16.1) 



with the initial condition N(x, t = 0) = 5(x), linx^j-oo N(x, t) = 0, where n is a diffusion constant and S(x) 
is the Dirac delta function. Then for the solution of (16.1) there holds the formula 



N(x,t) = —Hli 



1 

a\x\ 



(i,i),(i, !),(i, iy 



(i,i),(i,i),(i,i) 



(16.2) 



For details of the proof, the reader is referred to the original paper by Saxena, Mathai and Haubold (2006). 
Corollary 16.1. For the solution of the fraction reaction- diffusion equation 



^-N(x,t) = n ^BZNfat), 



(16.3) 
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with initial condition N(x, t = 0) = S(x) there holds the formula 



' \x\ 




_( v t)i 





(16.4) 



where a > 0. 



Note 16.1. It may be noted that (16.4) is a closed form representation of a Levy stable law. It is interesting 
to note that as a — > 2 the classical Gaussian solution is recovered since 



l ^ (_i)fc 

= (47r(r?t) « exp 



(i,i),(i,|)- 

(i,i),(i,|) . 
\x\ ■ 



2( V t)- 



2k+l 



(16.5) 
(16.6) 



17. Application to nonlinear waves 

It will be shown in this section that by the application of the inverse Laplace transforms of certain 
algebraic functions derived in Saxena, Mathai and Haubold (2006), we can establish the following theorem 
for nonlinear waves. 

Theorem 17.1. Consider the fractional reaction- diffusion equation 



D^N(x,t)+a Q D^N(x,t) 

= v 2 _ oa DlN{x, t) + ( 2 N(x, t) + 4>{x, t) 



(17.1) 



for x £ R, t > 0,0 < a < 1,0 < [3 < 1 with initial conditions N(x,0) — f(x), \im x ^± 00 N(x,t) — for 
x G R, where v 2 is a diffusion constant, ( is a constant which describes the nonlinearity in the system, and 
(f>(x, t) is nonlinear function which belongs to the area of reaction- diffusion, then there holds the following 
formula for the solution of (17.1). 



°° (—n) r f'°° 
^^Et 1 t< Q -«T(fc)exp(-fcE) 

X [*ES-0r+l(-&O +* a " /, ^S-«(r + l) + l(-W a )]d* 

" (-«r " 



E 

r=0 



27T 



^ ^a,a+(a— 0)r 



+ — 0)r— 1 

(-b( a )dk dC 



4>{k, t—Q exp(-ikx) 



(17.2) 



where a > (3 and Ep (■) is the generalized Mittag-Leffler function, defined by (11.1), and b = z^ 2 |/c| 7 — ( 2 . 



Proof. Applying the Laplace transform with respect to the time variable t and using the boundary condi- 
tions, we find that 



s a N(x, s) - s a - 1 f{x) + as fi N{x, s) - as^fix) 

= u 2 ^DlN(x, s) + ( 2 N{x, s) + f(x, s). 



(17.3) 
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If we apply the Fourier transform with respect to the space variable x to (17.3) it yields 



s a N*(k, s) - s a -\f*(k) + as f:S N*(k, s) - as^f^k) 

= -v 2 \kpN*(k, s) + ( 2 N*(k, s) + f*(k, s). (17.4) 

Solving for N* it gives 

N* (k s) _ (£!^j^)m±jlM fl75 ) 

^M)- s « + as0 + b (17 - 5) 

where b = f 2 |/c| 7 — Q 2 . For inverting (17.5) it is convenient to first invert the Laplace transform and then the 
Fourier transform. Inverting the Laplace transform with the help of the result Saxena et al. (2004a, eq(28)) 

-;t\ =t a - p Y(-a) r t ( - a - f3 '> r E r+ \ ( m +A-bt a ), (17.6) 

S a + as 13 + b J a.a+(a-(3)r-p+l\ H \ > 

' r=0 



where SR(a) > 0, SR(/3) > 0, SR(p) > 0, | ^ | < 1 and provided that the series in (17.6) is convergent, it yields 

oc- 

N*{k,t) = J2(- a Y' tia ~ f3)r f*( k ) 

r=0 

•x\E r+ l 1 * .A-W) +t a - f3 E r + l fl w . n . , (-6*°)] 

L a, (a— p)r+\ V / q,(q— p)(r+l) + l \ /J 

oo t 

+ £(-*) r / ^(fc,i-C)C a+(a - /3)r - 1 ^_ /3)r+a (-&C a )dC. (17.7) 

Finally, the inverse Fourier transform gives the desired solution (17.3). 

18. Generalized Mittag-Leffler type functions 

The multiindex (m-tuple) Mittag-Leffler function is defined in Kiryakova (2000) by means of the power 
series 

oo oo 

*u> W*) = E ^ = E rr r L + ky (is-i) 

k =o fe=o 1 1j=i 1 + Tj> 

Here to > 1 is an integer, pi,...,p m and pi,..., /* m are arbitrary real parameters. 

The following theorem is proved by Kiryakova (2000, p. 244) which shows that the multiindex Mittag- 
Leffler function is an entire function and also gives its asymptotic estimate, order and type. 

Theorem 18.1. For arbitrary sets of indices pi > 0, — oo < pi < oo,i = 1, ...,m the multiindex Mittag- 
Leffler function defined by (18.1) is an entire function of order 



m 1 11 1 
p = [V — I -1 , that is, - = K..H (18.2) 

and type 



i=1 Pi P Pi 



a= — )«...( — )pm. (18.3) 
P P 

Furthermore, for every positive e, the asymptotic estimate 

\E {1/piUm) (z)\<exp((a + e)\zn (18.4) 
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holds for \z\ > r (e),r (e) sufficiently large. 

It is interesting to note that for m = 2, (18.2) reduces to the generalized Mittag-LefHer function con- 
sidered by Dzhcrbashyan (1960) denoted by (fi pl , P2 (z; Hi, H2) and defined in the following form (Kiryakova, 
1994, Appendix) 



^(1/ P1 ,1/P2):(P1^2)( Z ) = ^Pl,Pl( Z \ Ml) M2) 



E 



> r(Mi + £)r( M2 + a 



(18.5) 



and shown to be an entire function of order 

P1P2 , , pi £2 p 2 

p = and type a = ( — ) pl ( — I 

P1+P2 P2' >1 

Relations between multiindex Mittag-Leffler function with H-function, generalized Wright function and other 
special functions are given by Kiryakova; for details, see the original papers Kiryakova (1999, 2000). Saxena, 
Kalla and Kiryakova (2003) investigated the relations between the multiindex Mittag-Leffler function and 
the Ricmann-Liouville fractional integrals and derivatives. The results derived are of general nature and give 
rise to a number of known as well as unknown results in the theory of generalized Mittag-Leffler functions, 
which serve as a backbone for the fractional calculus. Two interesting theorems established by Saxena et al. 
(2003) are described below. 

Theorem 18.2. Let a > 0, pi > 0, fn > 0,i = 1, ...,m and further, let Iff be the left-sided Riemann- 
Liouville fractional integral. Then there holds the relation 



(18.6) 



(V i ~ 1S (V^,(^ 1M )])(z) 



(l/Pi)>0*l+a.M2, — .Mm) 



(at Pi ). 



(18.7) 



Theorem 18.3. Let a > 0,pi > 0, /U, > 0, i = 1, ...,m and further, let I" be the right-sided Riemann- 
Liouville fractional integral. Then there holds the relation 



(I1[t-^- a E {1/piUlli) (at-n)])( x ) 

= x-^E (1/Pz) ^ 1+a ^ 2 ^.^ m) {at~ 



(18.8) 



Another generalization of the Mittag-Leffler function is recently given by Sharma (2008) in terms of the 
M-series defined by 



P Mq(ai, ...,a p ;h, ...,b q ;a;z) 



E 



(ai) r ...(a p ) r 



- (6i) r ...(6,) r r(ar + l) 



< 



1. 



(18.9) 



Remark 18.1. According to Saxena (2009), the M-series discussed by Sharma (2008) is not a new special 
function. It is, in disguise, a special case of the generalized Wright function p ^ q (z), which was introduced 
by Wright (1935), as shown below. 

24 



K p+ilp 



9+1 



(fli, 1), (a p , 1), (1, 1) _ 

(6i,l),...,(6„l),(l,a) 



E 



(a 1 ) r ...(a p ) r (l) r z r 
r ^(6i) r ...(6,) r r(ar + l) r! 

_ (fli) r ...(a p ) r r 

-2^( bl ) r ...( bq ) r T(ar + l) Z 

= p M^(ai,...,a p ;bi,...,b q ;a;z) 



(18.10) 



where 



k = 



rij=i r («i)' 



Fractional integration and fractional differentiation of the M-series are discussed by Sharma (2008). The two 
results proved in Sharma (2008) for the function defined by (18.9) are reproduced below. For 3?(i/) > 



( 7 o+[p M g( a i. -,(h;bu ...,b q ;a;z)])(x) 



Z 

r(i + z/) p+lM g+^ ai ' '"' ap ' 1 + v ' a ' z ) 



and for $t(v) < 



(D£ + |„M°(oi, a p ;h, ...,b q ;a; z)])(x) 



T(l^) P+ q+1 



19. Mittag-Leffler Distributions and Processes 

19.1. Mittag-Leffler Statistical Distribution and Its Properties 

A statistical distribution in terms of the Mittag-Leffler function E a {y) was defined by Pillai (1990) in 
terms of the distribution function or cumulative density function as follows: 

G v (y) = l-E a (-y a ) = ]T r(1 + ,0< <*<l,y>0 (19.1.1) 

fe=i 

and G y (y) = for y < 0. Differentiating on both sides with respect to y we obtain the density function f(y) 
as follows: 



f(y) = ^ G v(y) 



dy r(l + afc) 



E 



E 



fc +l„,l,„afc-l 



(-l) fc+1 afc?/ 
r(l + afc) 



r(l + afc) r(a + afc) 

= y Q_1 £a,a(-y a ), 0<a<l,y>0 
where E ai p(x) is the generalized Mittag-Leffler function. 



by replacing fc by k + 1 



(19.1.2) 



25 



It is straightforward to observe that for the density in (19.1.2) the distribution function is that in 
(19.1.1). The Laplace transform of the density in (19.1.2) is the following: 

/>oo /.oo 

L f (t)= e- tx f(x)dx= e'^x^E^-x^dx 
Jo Jo 

= (l + t a )-\\t a \ <1. (19.1.3) 

Note that (19.1.3) is a special case of the general class of Laplace transforms discussed in Section 2.3.7 
[Haubold and Mathai (2008)]. From (19.1.3) one can also note that there is a structural representation in 
terms of positive Levy distribution. A positive Levy random variable u > 0, with parameter a is such that 
the Laplace transform of the density of u > is given by e~* . That is, 

E[c- tu ] = e-*° (19.1.4) 

where E(-) denotes the expected value of (•) or the statistical expectation of (•). When a = 1 the random 
variable is degenerate with the density function 



, , _ J 1, for x = 1 
n\ x ) — | o, elsewhere. 



Consider an exponential random variable with density function 

= { relsewherl" < °° ^ <*> = ^ + ^ ^ 
and with the Laplace transform L f 1 (t) . 

Theorem 19.1.1 Let y > be a Levy random variable with Laplace transform as in (19.1.4) and let x 
and y be independently distributed. Then u — yx~ is distributed as a Mittag-Leffler random variable with 
Laplace transform as in (19.1.3). 

Proof. For establishing this result we will make use of a basic result on conditional expectations, which 
will be stated as a lemma. 

Lemma 19.1.1 For two random variables x and y having a joint distribution, 

E(x) = E[E(x\y)} (19.1.6) 

whenever all the expected values exist, where the inside expectation is taken in the conditional space of x 
given y and the outside expectation is taken in the marginal space of y. 

Now by applying (19.1.6) we have the following: Let the density of u be denoted by g(u). Then the 
Laplace transform of g is given by 

E[ C - < - txi ^\x}=e- tax - (19.1.8) 

But the right side of (19.1.8) is in the form of a Laplace transform of the density of x with parameter t a . 
Hence the expected value of the right side is 

L g (t) = (l+t a )- 1 (19.1.9) 

which establishes the result. From (19.1.8) one property is obvious. Suppose that we consider an arbitrary 
random variable y with the Laplace transform of the form 

L g (t) = e- [ * (t)1 (19.1.10) 
whenever the expected value exists, where 4>{t) be such that 

4>{tx») =x<f)(t),lim<f>(t) =0. 
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Then from (19.1.8) we have 

E[ e -( tx "^\x] = e - x ^W]. (19.1.11) 

Now, let x be an arbitrary positive random variable having Laplace transform, denoted by L x (t) where 
L x (t) = ip(t). Then from (19.1.9) we have 

L fl (t) = V#(t)]. (19.1.12) 

For example, if y is a random variable whose density has the Laplace transform, denoted by L y (t) = <j>(t), 
with 4>(tx~) — x<f>(i), and if x is a real random variable having the gamma density, 

fx{x) = X l ^ , x > 0, > 0,6 > (19.1.13) 

and f x (x) = elsewhere, and if x and y are statistically independently distributed and if u = yxi then the 
Laplace transform of the density of u, denoted by L u (t) is given by 

L u (t) = [l + 6{<P(t)}}-?. (19.1.14) 



Note 19.1.1. Since we did not put any restriction on the nature of the random variables, except that the 
expected values exist, the result in (19.1.12) holds whether the variables are continuous, discrete or mixed. 

Note 19.1.2. Observe that for the result in (19.1.12) to hold we need only the conditional Laplace transform 
of y given x be of the form in (19.1.11) and the marginal Laplace transform of x be ip(t). Then the result 
in (19.1.12) will hold. Thus statistical independence of x and y is not a basic requirement for the result in 
(19.1.12) to hold 

Thus from (19.1.13) we may write a particular case as 

z = yxi (19.1.15) 

where x is distributed as in (19.1.5) and y as in (19.1.4) then z will be distributed as in (19.1.3) or (19.1.9) 
when x and y are assumed to be independently distributed. 

Note 19.1.3 The representation of the Mittag-Lefner variable as well as the properties described on page 
1432 of Jayakumar (2003) and on page 53 of Jayakumar and Surcsh (2003) are to be rewritten and corrected 
because the exponential variable and Levy variable seem to be interchanged there. 

By taking the natural logarithms on both sides of (19.1.15) we have 

-lnx + lny = lnz. (19.1.16) 

a 

Then the first moment of lnz is available from (19.1.16) by computing £?[lnx] and E[lny}. But i£[lnx] is 
available from the following procedure: 

E[c- tlnx ] = E[c lnx "} = E[x- 1 } = / x^e-'dx 

Jo 

= T(l - t) for 3?(1 - t) > (19.1.17) 
which will be r(/3 - t)/T((3) for the density in (19.1.13). Hence 



E[lnx] = -±E[e-^*]\ t=0 = - - i)| t=0 . 



But 
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A r (i - 1) = -r(i - 1)^ inr(i - 1) = -r(i - t)V(i - *) (w.i.is) 

where tp(-) is the psi function of (•), see Mathai (1993) for details. Hence by taking the limits t — > 

Spna;] = -r(l - i)^(l - t)| t=0 = = 7 (19.1.19) 

where 7 is Euler's constant, see Mathai (1993) for details. 

19.2. Mellin-Barnes representation of the Mittag-Leffler density 

Consider the density function in (19.1.2). After writing in series form and then looking at the corre- 
sponding Mellin-Barnes representation we have the following: 

°° (— '\\ k T a-l+ak 

9{X ) = X ^E aA - X °) = gr ( l + k)^- f( - T - j - ) (19.2.1) 

1 r c+l °° 1 r(i - ^)r(i — — + — ) 

= — / ^ r-s «^a;- s ds,l - a < c < 1 (19.2.2 

2 « ,/c-ioo « r (! - s) 

[by expanding as the sum of residues at the poles of L(l — — + — )] 

= lv ^ ' ) x as - 1 d s = g(x),0<c 1 <l,0<a<l (19.2.3) 



by putting — — — = s\. Here the point s — is removable. By taking the Laplace transform of g(x) from 
(19.2.1) we have" 



a+ak-l -tx 



(—1) f° 

00 

= ^(-l) fe i-"- Qfe = t- Q (l +^")- 1 = (1 +t a )-\ \t a \ < 1. (19.2.4) 



fc=0 



19.2.1. Generalized Mittag-Leffler density 

Consider the generalized Mittag-Leffler function 

01 to - J_ V (-W-r + ^ -i-H* 

T( 7 ) ^ A;!r(afc + a 7 ) 
= a; a T- 1 £2 7 ,a(-» a ), « > 0, 7 > 0. (19.2.6) 
Laplace transform of gi (x) is the following: 



x m _ y (-i) fc (7)* r 



^cry+afc-lg-tx^ 



= V(-l) fe ^V a7 - afe = (1 + i")- 7 , < 1. (19.2.7) 

In fact, this is a special case of the general class of Laplace transforms connected with Mittag-Leffler function 
considered in Mathai, et al. (2006). 
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19.3. Mittag-Leffler Density as an H-function 

gi{x) of (19.2.6) can be written as a Mellin-Barnes integral and then as an H-function. 



si 0*0 



i i 

f(jj) 2^1 



± rrl,l 

m 



c+ioo 



T(s)T(r]-s) 



-ioo 


T(ar] — as) 






x a 






(0,l),(l-a?;,a)_ 



: °'J- 1 (a; 01 )-»ds,R(»y) > 0,0 < c < 3^(77) 



1 1 r 

~ aT( V ) 2m J c 
[by taking ar\ — 1 — as = — si] 

1 



ci+ioo 



r(« - 1 + ±)r(± - 
r(i- s ) 



ar(r?) 



171,1 

-"1.2 



(»?-i.i).(o,i) 



(19.3.1) 
(19.3.2) 

(19.3.3) 



Since g and g\ are represented as inverse Mellin transforms, in the Mellin-Barnes representation, one can 
obtain the (s — l)-th moments of g and g\ from (19.3.2). That is, 



M gi (s) = Eix*- 1 ) in 9l 



1 nv-i + mi 



r(7) 



aT(l - s) 



for 1 - a < SR(s) < 1, < a < 1, 77 > 0. 



M fl (t) ^Eix*- 1 ) in .9 

r(l-I + i)r(i-i) 

_ _^ a q / ^ q a ^ 

nli I - s) 

for 1 — a < jft(s) < 1, < a < 1, obtained by putting 77 = 1 in (19.3.4) also. Since 



r(i 

lim 



;) 



a-l r(i - s) 



= 1 



for a — > 1, (19.3.4) reduces to 



M gi (t) 



r(?7-l + s) fora^ 1. 



Its inverse Mellin transform is then 

gi = 



r(ri- l + s)a;- s ds 



1 



c r >- 1 e- x ,x> 0,r/ > 



(19.3.4) 



(19.3.5) 



(19.3.6) 



(19.3.7) 



T(t]) 2m J c _ ioo T(rj) 

which is the one-parameter gamma density and for r\ = 1 it reduces to the exponential density. Hence 
the generalized Mittag-Leffler density gi can be taken as an extension of a gamma density such as the 
one in (19.3.7) and the Mittag-Leffler density g as an extension of the exponential density for r\ = 1. Is 
there a structural representation for the random variable giving rise to the Laplace transform in (19.2.4) 
corresponding to (19.1.10)? The answer is in the affirmative and it is illustrated in (19.1.14). 

Note 19.3.1. Pillai (1990, Theorem 2.2), Lin (1998, Lemma 3) and others list the p-th moment of the 
Mittag-Leffler density g in (19.1.2) as follows: 



r(i - £)r(i + £) 
E(x") = r(i- P ) ' ~ a K <a<1 - 

Therefore 
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i r(i - i)rri - i + -s-) 

V g g / V a a / 

(1 - s)r(i - s) a r(i - s) 

which is the expression in (19.3.5). Hence the two expressions are one and the same. 

Note 19.3.2. If y = ax, a > and if x has a Mittag-Lefner distribution then the density of y can also be 
represented as a Mittag-Lefner function with the Laplace transform 

L x (t) = (l + t*)- 1 =>L y (t) = (l + (at) a )-\a>0,\(at) a \ < 1. 
Note 19.3.3. From the representation that 

E(x h ) = r(iA) ' ~" < m) <a<1 

we have 

Further, g(x) in (19.2.1) is a non-negative function for all x, with 

E(x h ) = / x h g{x)dx = 1 for h = 0. 
Jo 

Hence 5(2;) is a density function for a positive random variable x. Note that from the series form for the 
Mittag-LefHer function it is not possible to show that J °° g{x)dx = 1 directly. 

19.4. Structural Representation of the Generalized Mittag-Lefner Variable 

Let u be the random variable corresponding to the Laplace transform (19.2.7) with t a replaced by 5t a 
and 7 by 77. Let u be a positive Levy variable with the Laplace transform e _t , < a < 1 and let v be a 
gamma random variable with parameters 77 and 5 or with he Laplace transform (1 + St)~ v ,ri > 0, S > 0. Let 
u and v be statistically independently distributed. 

Lemma 19.4.1. Letu,v as defined above. Then 

w ~ uvi (19.4.1) 

where w is a generalized Mittag-Leffler variable with Laplace transform (l + St 01 )^ 1 ^ , \St a \ < 1, where <~ means 
'distributed as ' or both sides have the same distribution. 

Proof. Denoting the Laplace transform of the density of w by L w (t) and treating it as an expected value 

L w (t) = E[e~ tviu ] - E[E[e- (tviu) \v]} 

= E[c- ( - tvi ') a } = E[c- tav ] = (1 + St a y n (19.4.2) 

from the Laplace transform of a gamma variable. This establishes the result. 

From the structural representation in (19.4.1), taking the Mcllin transforms and writing as expected 
values, we have 

Eiw 3 - 1 ) = Efa'-^Eiv*)'- 1 (19.4.2) 
due to statistical independence of u and v. The left side is available from (19.3.4) as 
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E{x s ^) 



r(2 - S ) 
(- — -)r(- — -)r(i — - + -) 

x a g /v g g /v a a / 



v g g 



g a ' 



Y(n - - + -)V(- - iUV 

TTl„,,8— 1\ _ V' g ' g / V g g ^ /in i o\ 

^ } [fpR ' ( ) 

Let us compute _E[u~] s_1 from the gamma density. That is, 

£7[t,i]-i = ^ f^Y-^e-Uv = ^"^^ (19-4.4) 

for 3J(s) > 1 — 0:77, < a < l,rf > 0. Comparing (19.4.3) and (19.4.4) we have the (s — l)-th moment of a 
Levy variable 

r(- - -) rfi + - - -) 

E[US ~ 1] = ar(l -J) = T(2-s) MS) < ° K a ~ L (19 ' 4 - 5) 
Note 19.4.1. Lin (1998) gives the p-th moment of a Levy variable with parameter a as 

r(i - £-) 

E[u p ] = T(l _ a p y (19-4.6) 



Hence for p = s — 1 we have 



m + i-i) (i_«)r(i-i) 

L(2 - s) " ~ (1 - s)T(l - s) 
T(I - ») 

This is (19.4.5) and hence both the representations are one and the same. 
Hence the Levy density, denoted by g2{u), can be written as 



i i-c+ico Y(— — —) 

•»(->-5sjLafr4"- <1 '- < 19A7) 



Its Laplace transform is then 



1 /-C+ZOO TV.!. _s \ /*oo 

92 W 2«i c _ 40O ar(l-s)7o 

= -L / ±T(- - -)t- 1+s ds 
2 ™ Jc-ioo a a a 



c— too 
' Jci— too " 

by making the substitution — 1+s = — si. Then evaluating as the sum of the residues at ^ = —u, v = 0, 1, 2, ... 
we have 

L ^W = ES/ t a " = e-* a . (19.4.8) 

This verifies the result about the Laplace transform of the positive Levy variable with parameter a. Note 
that when a = 1, (19.4.8) gives the Laplace transform of a degenerate random variable taking the value 1 
with probability 1. 

Mcllin convolution of certain Levy variables can be seen to be again a Levy variable. 
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Lemma 19.4.2. Let Xj be a positive Levy variable with parameter a>j,0 < aj < 1 and let x\,...,x p be 
statistically independently distributed. Then 



is distributed as a Levy variable with parameter aia 2 --.a p . 
Proof. From (19.4.5) 



E[e~ tx i] 


= e~ tj 


E[e~ tu } 


= E[e~ 




= E\e- 



(III" 1 ...I° 1 '"° p X ]_ fpf fpl a -tx 1 ...X p 1 " 



= E[E[e -t Xl ... Xp \ X2 _ Xp \] 



Repeated application of the conditional argument gives the final result as 

E[e~ tu ) = e~* , a — a\a2.--a pi < a\...a p < 1 

which means that u is distributed as a Levy with parameter ai...a p . 

From the representation in (19.4.1) we can compute the moments of the natural logarithms of Mittag- 
Lcfncr, Levy and gamma variables. 

w = ra« =^> In w = lnu + — In v. (19.4.9) 

a 

But from (19.4.3) and (19.4.6) we have the h-th moments of u and v given by 

E[u h ] = r^I^ .KW < « < 1 (19.4.10) 



and 



From (19.4.3) 



E[vi] h = Efo±llfl , ftfo + > 0. (19.4.11) 



But for a positive random variable z 

E[z h ] =E[c lDzh }. 

Hence 

^E[z h ]\ h=0 = E[\nzc hl * z ] h=0 = E[\nz}. 
Therefore from (19.4.9) to (19.4.12) we have the following: 

d f«-r(„ + £)r(i-£)' 



E[lnw] 



dh T{n)T(l - h) 



\h=0 



= - -^(1)+^(1) + -lnJ (19.4.13) 

a a a 
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by taking the logarithmic derivative, where tp is a psi function, see, for example Mathai (1993). 



Efrvi] = A {^^} U=o = >) + iln* (19.4.14) 



£[lni;] = V(?7) + fn<5 

and 

^ = A{^|}l^o = ->) + ^(l) (U.4.15) 
where V(l) = ^7 where 7 is the Euler's constant. 

Note 19.4.2. The relations on the expected values of the logarithms of Mittag-Leffler variable, positive 
Levy variable and exponential variable, given on page 1432 of Jayakumar (2003), where i) = 1, arc not 
correct. 

19.5. A Pathway from Mittag-Leffler Distribution to Positive Levy Distribution 

Consider the function 

^) = E mW 1 n ^ + , ,r?>0,a>0,0<a<l 
^ k\T{ari + ak) (a~) ar i+ ak 

x cm)-1 x <y 



E r > ( ) 



Thus x — a <=,y where y is a generalized Mittag-Leffler variable. The Laplace transform of / is given by the 
following: 



fc=0 



k\a r >+ k J Q T(ar] + ak) 



k=0 

If r] is replaced by and a by a(q — 1) with q > 1 then we have a Laplace transform 

L f (t) = [1 + a(q-l)t a }-^,q> 1 (19.5.2) 

If q — > 1 + then 

L / (i)^e-<=L A (t) (19.5.3) 

which is the Laplace transform of a constant multiple of a positive Levy variable with parameter a. Thus q 
here creates a pathway of going from the general Mittag-Leffler density / to a positive Levy density /1 with 
parameter a, the multiplying constant being (077)". For a discussion of a general rectangular matrix- variate 
pathway model see Mathai (2005). The result in (19.5.3) can be put in a more general setting. Consider an 
arbitrary real random variable y with the Laplace transform, denoted by L y (t), and given by 

L v (t) = e-*® (19.5.4) 
where 4>{t) is a function such that 0(te 7 ) = x<fi(t), <f>(t) > 0, lim^o 4>{t) = for some real positive 7. Let 
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u = yx 



(19.5.5) 



where x and y are independently distributed with y having the Laplace transform in (19.5.4) and x having a 
two-parameter gamma density with shape parameter (3 and scale parameter <5 or with the Laplace transform 

L x (t) = (l + St)- . (19.5.6) 
Now consider the Laplace transform of u in (19.5.5), denoted by L u (t). Then 

L u (t) = E[c- tu ] = E[e- tyx "] = E[E[e tyx "\x\] 
= ^[e-^^ 7 )] = E[e- X ^] 

from the assumed property of 4>{t) 

If S is replaced by 5(q — 1) and (3 by with q > 1 then we get a path through q 

L u (t) = [1 + S(q - l)4>{t)]-& - e"WW = e^^l 
If 4>(t) = t a ,0 < a < 1 then 

which means that u goes to a constant multiple of a positive Levy variable with parameter a, the constant 
being (<5/3)t. 

19.6. Linnik or a-Laplace Distribution 

A Linnik random variable is defined as that real scalar random variable whose characteristic function 
is given by 

4>{t) = - — j-r^,0 < a < 2,-oo < t < oo. (19.6.1) 

For a = 2, (19.6.1) corresponds to the characteristic function of a Laplace random variable and hence 
Pillai (1995) called the distribution corresponding to (19.6.1) as the a-Laplace distribution. For positive 
variable, (19.6.1) reduces to the characteristic function of a Mittag-Lefner variable. Infinite divisibility, 
characterizations, other properties and related materials may be seen from the review paper Jayakumar and 
Suresh (2003) and the many references therein, Pakes (1998) and Mainardi and Pagnini (2008). Multivariate 
generalization of Mittag-Leffler and Linnik distributions may be seen from Lim and Teo (2009). Since the 
steps for deriving results on Linnik distribution are parallel to those of the Mittag-Leffler variable, further 
discussion of Linnik distribution is omitted. 

19.7. Multivariable Generalization of Mittag-Leffler, Linnik and Levy Distributions 

A multivariate Linnik distribution can be defined in terms of a multivariate Levy vector. Let T" = 
(ti, ...,t p ),X' = (xi, ...,x p ), prime denoting the transpose. A vector variable having positive Levy distribu- 
tion is given by the characteristic function 

E[e iT ' x ] = c -( T ' ST ) f , < a < 2, (19.7.1) 
where E = E' > is a real positive definite p x p matrix. Consider the representation 

u = y°X (19.7.2) 

where the p x 1 vector A, having a multivariable Levy distribution with parameter a, and y a real scalar 
gamma random variable with the parameters 8 and (3, are independently distributed. Then the characteristic 
function of the random vector variable u is given by the following: 



(19.5.7) 

. That is, when q — > 1 + , 
(19.5.8) 
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£[e ,s ° T ' x ] = E[E[e iVaT ' x }\y] 

= E[e-y^ T '^ % } = [1 + 5\T'Y,T\%]~P . (19.7.3) 

Then the distribution of u, with the characteristic function in (19.7.3) is called a vector-variable Linnik 
distribution. Some properties of this distribution are given in Lim and Teo (2009). 

20. Mittag-Leffler Stochastic Processes 

The stochastic process {x(t),t > 0} having stationary independent increment with x(0) = and x(l) 
having the Laplace transform 

L x(1) (A) = (1 + A a )"\0<a< 1,A>0, 

which is the Laplace transform of a Mittag-Leffler random variable, is called the Mittag-Leffler stochastic 
process. Then the Laplace transform of x(t), denoted by L X ^(X), is given by 

L x{t) (X) = [(1 + A")- 1 ]* = [1 + A*]"*. (20.1.1) 

The density corresponding to the Laplace transform (20.1.1) or the density of x(t) is then available as the 
following: 



Oik-\-at— 1 



k\ T(ak + at) 

= x^^El , at (-x a ), < a < 1, x > 0, t > 0. (20.1.2) 
The distribution function of x(t) is given by 



t)T(ak + at) 

-±t2. r <f\* „ n< a < M >0. (20.1.3) 

f-^ fc! rm r(l + ak + at) ' - ' v ' 



This form is given by Pillai (1990) and by his students. 
20.2. Linear First Order Autoregressive processes 

Consider the stochastic process 



e„ , with probability p, < p < 1 

e„ + as„_i with probability 1— p, < a < 1. 



(20.2.1) 



Let the sequence {e„} be independently and identically distributed with Laplace transform L e (X) and let 
{x n } be identically distributed with Laplace transform L X (X). From the representation in (20.2.1) 



L Xn (X) = P L e (X) + (1 - p)L e (A)L Xn _ 1 (aA) 

Therefore 



el J p+a-pJL^^aA) p+(l-p)L x (a\) [ U ' ' J 
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assuming stationarity. When p = 0, 



ie(A)-i^) ,0<a<l (20.2.3) 

which defines class L distributions, for all a, < a < 1. When p = 0, (20.2.3) implies that the innovation 
sequence {e„} belongs to class L distributions. Then (20.2.3) can lead to two autoregressive situations, 
the first order exponential autoregressive process EAR(1) and the first order Mittag-LefHer autoregressive 
process MLAR(l). 

Concluding remarks 

The various Mittag-LefHer functions discussed in this paper will be useful for investigators in various 
disciplines of applied sciences and engineering. The importance of Mittag-LefHer function in physics is 
steadily increasing, ft is simply said that deviations of physical phenomena from exponential behavior could 
be governed by physical laws through Mittag-LefHer functions (power- law). Currently more and more such 
phenomena are discovered and studied. 

It is particularly important for the disciplines of stochastic systems, dynamical systems theory, and 
disordered systems. Eventually, it is believed that all these new research results will lead to the discovery of 
truly nonequilibrium statistical mechanics. This is statistical mechanics beyond Boltzmann and Gibbs. This 
nonequilibrium statistical mechanics will focus on entropy production, reaction, difiusion, reaction-diffusion, 
etc and may be governed by fractional calculus. 

Right now, fractional calculus and H-function (Mittag-LefHer function) are very important in research 
in physics. 
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